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^^ ' In an abelian gauge theory, the Coulonibic potential between two static 

f^ ■ charges is obtained most directly when a correct separation between gauge- 

r^ . invariant and gauge degrees of freedom is made. This motivates a similar 

separation in the nonabelian theory. When a careful identification of the 

Hilbert space is made, along with the proper analyticity requirements, it is 

Q^' then possible to find the appropriate wavefunctionals describing heavy color 

^ ■ sources in the theory. This treatment is consistent with, and realizes in a 

simple way the center symmetry Z^ of SU{N) gauge theories. 
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I. INTRODUCTION 

Anyone who has given a certain amount of thought to different treatments of quantum 
electrodynamics comes inevitably to the conclusion that the reason it is a rather straight- 
forward theory to work out is, in one guise or another, always the same: at some point, one 
needs to clearly separate what are physical degrees of freedom from what are gauge degrees 
of freedom, and this separation is absolutely clear in an abelian gauge theory. Unfortunately, 
no such luck persists with nonabelian gauge theories. The nonlinear nature of the theory 
obscures considerably this separation of physical from gauge components. One might never- 
theless wonder whether a similar systematic separation of degrees of freedom is possible in 
the nonabelian theory as well and, if so, whether it may shed some light on different aspects 
of the theory in any way. 

In previous work |jl[, we have given a description of the Hilbert space of Yang-Mills 
theory in terms of gauge invariant variables which, as we show in what follows, does achieve 
this objective without tremendous technical effort. The specific issue we address here is 
the form state functionals of the theory should take in the presence of two static color 
sources. Firstly, it is of course the nonzero charge sectors that exploit this separation of 
degrees of freedom. Moreover, through the static charge-anti-charge state one may also 
study the confining potentials generated by the theory. The central points guiding the 
search for these wavefunctionals are the above mentioned separation of degrees of freedom, 
and a proper identification of the Hilbert space of the theory. We take the view here that 
the only two requirements we are allowed to make are that wavefunctionals satisfy Gauss' 
law and that they be analytic functions in the Hilbert space of the theory. We find a 
posteriori that wavefunctionals for two static color sources satisfying these two requirements 
necessarily contain a dependence on some curve connecting the two charges if they belong to 
a fundamental representation, a behavior that is markedly different from both the abelian 
case and the case of nonabelian adjoint sources. The reason this behavior comes about 
lies in the correct identification of the Hilbert space in the nonabelian theory, which in 
turn comes from the observation that a pure SU{N) gauge theory is in reality only an 
SU{N)/Z]\r gauge theory because of the well-known fact that the gluons, which belong to 
the adjoint representation, are insensitive to the center transformations, Z^. Specifically for 
SU{2)/Z2 = 5*0(3), the description of the Hilbert space becomes particularly clear, through 
the use of Euler angles and variables which simply rotate under gauge transformations, 
and this is the case we will be concerned with here, although the extension to SU{N > 2) 
should not present any conceptual difficulties. Analyticity and periodicity requirements 
in the 5*0(3) domain will naturally lead to a simple and explicit realization of the center 
symmetry, through 't Hooft's operator algebra ^. At the same time, Wilson loops are not 
gauge invariant under transformations generating vortices through them, and this fact would 
be entirely missed if the gauge group were considered to be SU{2) rather than 5*0(3), unless 
the nontrivial Z2 behavior were inserted by hand, e.g., through boundary conditions 0. 

We begin by treating the abelian case, since this will not only motivate the correct sepa- 
ration of degrees of freedom in the nonabelian theory, but will also indicate some cautionary 
remarks concerning the physical interpretation commonly ascribed to Wilson loops, both in 
the abelian and the nonabelian cases. The standard procedure to obtain the potential for 



two static charges in the abehan theory is to calculate the vacuum expectation value of the 
Wilson loop 0] 

However, extracting a Coulomb law from Eq. (^ requires some nontrivial steps. The integral 
contains divergences, and thus if one computes it for, say, a rectangular loop of sides L and 
T, the exact result is a function of AL and AT, symmetric under L ^^ T interchange, where 
A is a momentum cutoff. The double integral then yields: 

ciA(L + T) + C2(L/T + T/L) + ... (2) 

(the exact result will be presented below). The Coulomb law is buried in Eqs. (|l],H) accord- 
ing to well-known arguments: one must first interpret the expectation value to mean the 
amplitude for creation of a charge-anti-charge pair, separated by a distance L, out of the 
vacuum during some finite time T. Then, one must argue away the first term as a self-energy 
contribution, and let T become large in order to isolate the behavior exp{—iT/L). Some 
further argumentation then allows one to identify the 1/L in this latter expression as the po- 
tential between the charges running in the loop. Altogether, we believe this entire procedure 
to be a bit too much calculation and argumentation. In a Hamiltonian formalism, the same 
result can be seen to follow in a much simpler fashion. We will also show, however, that 
unless one makes a judicious separation of gauge and gauge-invariant variables, one may 
well erroneously "derive" confinement in this theory. The lessons learned in this exercise 
will then serve as a guide in treating the nonabelian theory. In particular, we will show that 
the form wavefunctionals should take is not the one directly connected to Wilson loops, and 
sometimes indicated in the literature. 



II. (NO) CONFINEMENT IN THE ABELIAN THEORY 

Here and in what follows, we consider the Hamiltonian quantization of gauge theories in 
"temporal" gauge (in the abelian theory: Aq = 0), and in the Schrodinger representation. 
States of the theory are wavefunctionals ^[A] which are eigenf unctions of the Hamiltonian, 
with the further requirement that they satisfy Gauss' law, which in the abelian theory reads 

d ^*[^] -0 (3) 

If two static charges, of opposite unit charge, for definiteness, are placed at points xq 
and Xi, states of the theory will be modified, and will in particular solve a modified Gauss 
law: 

dk 1^ = -i{6^'\x - xi) - 5(3)(x - xo))nA] . (4) 

If a functional $[^4] satisfies Eq. (^, then the functional 

^[A] =e'C''"'"^'"^"^$[A] (5) 



will satisfy Eq. (^). In the above, the line integral is taken along some curve C connecting 
Xq and Xi, and for simplicity we do not indicate explicitly the dependence of \1/ on xq, Xi or 
C. The resemblance between this functional and the Wilson loop operator does suggest the 
above to be the natural candidate for the state with opposite charges at Xq and Xi, if $ is 
the vacuum state. The electric field expectation on ^ is obtained from 

^(^)*14 - \ j^,m - 1: ''y'^'^U. - y) m - .'/;■ ™|igL . (c, 

If we assume $ is real - we might take, for instance, $ to be indeed the vacuum state - 
there is no cross term in squaring this expression and we are led to the following electric 
energy expectation: 

j d?x<E\x)>^ = j d\<E\x)>^+ r dy"^ H dz'^d^^^y - z) (7) 

J J Jxq J Xq 

d^x<E\x)>^+2L6^^\0) , 



where L is the length of the curve C, and ^'^^•'(O) is the inverse thickness of the curve, to 
be made finite by some suitable renormalization. This is minimized for C a straight line 
between xq and Xi, and shows an electric energy growing linearly with the separation of the 
charges. In other words, confinement of charge as the term is usually understood. 

What have we done wrong? We have simply chosen a very bad functional: taking a real 
<l> has eliminated the cross term in the (squared) electric energy, and has given a highly 
divergent contribution. It is not that the derivation is wrong, but there is a much better 
wavefunctional around, leading to the exact solution, and a smaller electric energy. We find 
it by separating the gauge variables from the gauge invariant variables: 

Ai{x) = Aj{x) + diA{x) , (8) 

with diAj{x) = 0, and 

A(.) = --L/,5,^. (9) 

An J \x — y\ 

Inserted into Eq. (^, it leads to 

^[A] = e^(^(^i)-^(^o)) [e'^-'o' '^^™^™^^^$[A]] (10) 

^ gi(A(a;i)-A(a;o))|,Ml 

where ^[A] is gauge invariant and thus still solves the source- free Gauss law, Eq. (^). The 
electric field now reads 

1 -^mA] = (^-^^^ - ^-^^^] mA] - ,e^(A(^'i)-A(-o))iM (11) 

= d, (—1 _ —1 -] ^[A] - ,e^(A(-0-A(.o))i^ . 

\47r|x — Xi| 47r|x — Xoly oAi[x) 



In squaring this expression to obtain the electric energy, the cross term vanishes without 
any further assumptions since $ is gauge invariant and thus satisfies Eq. (j^). The electric 
energy then reads: 

dh <E'^(x)><i,= I d\<E'^ix)>^-2 — ■ r + s.e. , (12) 

J 47r|a;o — xil 

where "s.e." means the self-energy terms for both sources. Now we see the correct result, 
no confinement, just the Coulomb interaction. Whereas previously our choice of prefactor 
in Eq. (|]) led to the singular quantity 5 Ai{x) / 5 Aj{y) oc 5^'^\x — y), the second choice led 
on the other hand to the smoother quantity 6A{x)/6Ai{x). A seemingly "natural" choice 
of $ {i.e., $ the vacuum wavefunctional) then did not allow for the cancellation of this 
extraneous divergence. 

We thus learn that when we introduce a prefactor into wavefunctionals in order to solve 
the Gauss law with sources, the presence of the full vector potential in the line integral will 
lead to an artificially divergent quantity; one should rather attempt to introduce only the 
gauge part of the vector potential. In the abelian case, the gauge part of the vector potential 
appears in such a way as to effectively "disconnect" the two end points where the charges 
are located, without introducing any discontinuities (cf. Eq. (|13)); we may suspect that 
in a correct treatment of the nonabelian theory with fundamental sources this should not 
happen, as in that case one does expect a "line" to be present connecting the two sources. 

It is possible in fact to derive the correct result in a simpler way by directly solving the 
local form of Gauss' law, because it makes no premature use of a path: using the appropriate 
separation of gauge and physical variables, Eq. (§), it follows that 

^' JM^) ^ ~5A(^ ' ^^^^ 

which in turn leads to, using Eq. (|), 

-^ = t{6^'\x - xi) - 6^'\x - xo))^ , (14) 

dA{x) 

which immediately gives the correct solution: 

Further minimization in the gauge invariant sector tells us that $o is the vacuum wavefunc- 
tional. Incidentally, we note that this is not at all a new result: it can in fact be traced back 
to Dirac 0. 

The separation of gauge and gauge invariant variables has clearly indicated that it is 
the longitudinal degrees of freedom which are responsible for the potential between static 
charges, and once that separation is made. Gauss' law by itself allows for the computation 
of the potential. With that in mind, one understands why the same result is more awkward 
to extract from the Wilson loop, since that is an expectation value of a gauge invariant 
operator in the vacuum sector of the theory, which only depends on transverse variables. 
We have furthermore learned that the wavefunctional which is suggested by the Wilson 
loop, and which might seem at first more natural (Eq. (^) in fact is incorrect, essentially 



because it introduces all of the gauge field degrees of freedom, transverse and longitudinal, 
propagating along the line connecting the two charges. Given this result, however, one might 
now wonder: what has conspired to allow the Wilson loop expectation to also indicate the 
correct result, albeit with more effort? In order to answer this question, we first observe 
that the wavefunctional in Eq. (O) is an exact eigenfunction of the Hamiltonian operator 
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(16) 



with eigenvalue foA''\^ + wiA + t;2/|xo — Xi\, where Vi are constants, A is a momentum cutoff, 
V is the volume of space, the first term is a vacuum zero-point energy, the second term 
is a self-energy, and the last term is the Coulomb energy. Meanwhile, it is also simple to 
verify explicitly that the Wilson line, Eq. (^ (with $ the vacuum functional), is not an 
eigenstate of H. In order to make contact between the Wilson loop and our calculations, we 
may choose to calculate the former in the gauge Aq = 0, for a loop tracing a rectangle, with 
two fixed-time sides of length L and two sides of length T along the time direction. Due 
to the gauge choice, the expectation value only receives contributions from the spacelike 
sides. Furthermore, each of these sides corresponds to a state given by a Wilson line acting 
on the vacuum state, Eq. (^. Since one is at time T and the other at time 0, we may 
insert the operator exp —iTH in the middle in order to bring these states to the same time, 
and this operator can be further decomposed into its matrix elements on energy eigenstates. 
Then, the lowest nonvanishing overlap will give the leading exponential behavior of the loop. 
Because the Wilson line belongs to the sector with opposite charges at xq and xi, the lowest 
Hamiltonian eigenstate with a nonvanishing overlap will be precisely the eigenstate given 



by Eq. (pTSP, so that the correct Coulomb behavior is obtained. If the Wilson line only had 
components in the higher modes within this two-charge sector, the expectation value of the 
Wilson loop would be even further removed from the Coulomb law. 

In order to explicitly confirm the relevance of this overlap, we have performed the exact 
computation of the Wilson loop expectation value for a rectangle of sides Li and L2 along 
two spatial directions. It is a lengthy computation with, however, a rather compact result: 
we find 
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where sq — ^ is a cutoff parameter (of dimension L^) used to regulate the propagator in a 
Schwinger parametrization, and 7 is the Euler constant (we have also ommitted two terms 
vanishing in the limit Sq — ^ 0). To orient the loop along the time direction and extract the 



Coulomb law, we take Li >> L2 and then Li 
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-iT, and are left with 
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(18) 



The first term is the self-energy contribution, the second one is the Coulomb law, and the 
last term is the (modulus squared) overlap of the Wilson line with the correct eigenstate. 



Eq. ([T5|). As L2 is made larger, this overlap becomes smaller, which is expected, as for larger 
separations the Wilson line becomes more and more energetic, and thus further removed from 
the true charge-anti-charge state. Furthermore, we also see that, as the cutoff is removed, this 
overlap becomes vanishing. One may understand this as a renormalization of the composite 
operator corresponding to the Wilson line, but however one wants to interpret it, it does 
mean that the overlap does indeed vanish in the limit sq —>■ 0. 

As long as the Wilson line has a nonzero overlap with the correct wavefunctional for two 
charges, the Wilson loop expectation value should in the end lead to the correct potential 
between static charges. If, for instance, on a lattice calculation one is not forced to take the 
analogous limit Sq ^ in order to extract the appropriate scaling behaviors, then such a 
nonzero overlap may well be achieved. At the same time, however, we have seen that there 
are caveats to the full computation, even in the simpler abelian case where an exact result 
can be found, and in general Wilson loops should not be taken as any direct indication of 
what the correct wavefunctional is for a pair of static charges. Some of the literature on this 
subject may be a bit misleading, as sometimes it implicitly appears to assume the analogous 
form of Eq. (|^) for the wavefunctional in the nonabelian case, and in some instances such a 
wavefunctional is explicitly written down 0]. 

III. THE NONABELIAN THEORY 

Much of what we do here will be valid for general SU{N) gauge group, but in order to 
keep the discussion sufficiently simple and explicit, we will treat the SU{2) theory. In the 
presence of two static color sources at xq and xi, state wavefunctionals take the form 

^ap[A] , (19) 

where a and (3 are indices in the representation of the sources at Xi and Xq, respectively. 
Under finite gauge transformations, these states transform as 

^^^[A] -^ Qu'^aM] = '^apl^'A] = Uo.a'{x,)U;p,{xo)^^>p,[A] . (20) 

For notational simplicity, we have written the same "[/" everywhere, but of course it should 
be kept in mind that A transforms in the adjoint, while the f/'s at Xq and Xi belong to 
whatever representation the sources are in. Here we are mostly interested in the case in 
which both sources are either in the fundamental or in the adjoint. A physical wavefunctional 
should transform as the above under gauge transformations, and the infinitesimal statement 
of such a transformation is Gauss' law in the presence of sources: 

(g'^ix) + i5^'^\x - xi)A« - iS^^\x - xo)A-)^^^^,^, vi>^,^, = , (21) 

where A" is an SU{2) generator, either in the fundamental or in the adjoint, and 

with Df'^ = diS"''' + e"''^^ A'^ , is the infinitesimal generator of gauge transformations on function- 
als of A (in the above, unit matrices in aa', /3/3'-space are also understood in the appropriate 
places) . 



As previously, one solution of this constraint is gotten by attaching a Wilson line prefactor 
to a gauge invariant functional $: 

^^p[A] = P fe^ro''^"'"^'"^"^') $[A] , (23) 

V /a/3 

where P stands for the usual path ordering prescription, and A^ix) = y4^(x)A" is a matrix 
in the algebra of SU{2) in the representation of the sources. This functional transforms as 
in Eq. (pOD, and consequently satisfies Gauss' law in the presence of sources. Furthermore, 
it is a single- valued, analytic function of the vector potential in the S0{3) domain of the 
variables, irrespective of whether A° is a single- valued representation or not. However, like 
in the abelian case, it leads to an extraneous divergence when calculating the electric field 
energy, due to the singular term 5Ai{x)/5Aj{y) oc 5^^\x — y), and we must discard it. 

What we should seek, rather, is the equivalent of the abelian term 6A{x)/6Ai{x) in the 
electric energy, while at the same time preserving the transformation property Eq. (^) and 
single- valuedness in the 5*0(3) domain of the variables. 

The abelian example suggests that we attempt to solve Gauss' law in local form, since 
it led most simply to the correct result in that case. In order to do so, we have found 
it most convenient to utilize a set of gauge covariant variables that we have described in 
previous work |jl|. These new variables, uf, replace the vector potential through the following 
t r ansf or mat ion : 

e'^'^Djul = e'^\djul + e^^^A^^ul) = . (24) 

Unlike the vector potential, which transforms inhomogeneously, these new variables simply 
rotate under gauge transformations: 

n^\ti^)nu = u''\x)u1ix) , (25) 

with U"'^ a 3 X 3 real orthogonal matrix. From these, furthermore, gauge invariant variables 
can easily be obtained by contraction in color indices: gij = u'^u'^. The nine degrees of 
freedom of A^ are thus clearly split between six gauge invariant degrees of freedom in g^j 
and three S0{3) (rather than SU{2)) gauge angles. For our considerations here, it will 
suffice to consider the variables u";. Taken as three vectors in color space, for i = 1,2,3, they 
uniquely define a tetrahedron at each point of IR^, with vertices at the origin and at the 
ends of the vectors. Variations in the gauge invariant variables change the shape or the size 
of the tetrahedron, while gauge variations simply rotate it. In terms of the new variables, 
the Gauss law generator reads 

To satisfy Eq. (^) in the simplest possible way, one may take a wavefunctional 

^„/3[<(a;);<(xo),<(xi)] , (27) 

which is a functional of u^{x) throughout space, but a regular function of m"(xo) and u^{xi). 
Then, 
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6uf{x) 6uf{x) 



+ 6^^\x - xo) + (5^^^(x - xi) 



duf{xo) 



duf{xi) 



(28) 



where the functional derivative on the r.h.s. is taken with 'u"(xo) and m"(xi) fixed. Gauss' 
law becomes 



at Xq^ 



e'^'X'la^o) 



e'""'u'(x^: 
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iK'^* ^«« = 
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at Xi, and 



5u°;{x) 
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(29) 



(30) 



(31) 



MO,Ml 



If we associate Euler angles (a,/3, 7)0 and (a,/?, 7)1 to each of the tetrahedra u°'{xq) and 
m"(xi), the rotation generators in Eqs. (PP|,PD|) become the angular momentum operators 



J, = -i 



d 



da 



(32) 



±— ^(^lA- 3- — 
dj3 sin /5 I i c?7 i 9a; 



at Xf) and xi, with J± = Ji ± zJ2- 



The solution to Eqs. (|29| , |30D is well-known from the representation theory of SU{2). If 
the A" belong to the spin-j representation, it is given by the Wigner T)^^^ functions, which 
satisfy: 



rv 



(i) 



mn nm' ' 



(33) 



where m, m', n span integer values from —j to j. In the case of interest for us, j = 1/2, A" 
are the Pauli matrices t"/2, and the complete solution to Eqs. ( ^9| - |3l| ) is given by: 



^«/^ = ^i'rH(«,/5,7)i)^{.f ^*((«,/5,7)o)$"^"[^d 



(34) 



where $™"[5'jj] are four arbitrary gauge-invariant functionals (and thus functionals of gij 
only), with no particular symmetry associated to the indices m,n. These indices m,n live 
in the same (spin-1/2) space as a,f3, however, they are denoted in this way precisely to 
indicate the lack of any symmetry covariance associated to them, as far as Gauss' law is 
concerned. The "D^^/^) matrices appearing above are given by: 



-e 2 



(«-7) si 



p(i/2)(«.^)=fe-t(-+-)cos/?/2 
"'^ ^ '^'^^ 1^ et("-^) sin/3/2 ei("+^) cos/5/2 



sin/?/2 



(35) 



One more time, it seems we have been able to solve Gauss' law in the presence of two 
sources. Furthermore, it also seems we have eliminated the problem of introducing "too 



8 



much" of the electric field propagating between the two points. In fact, in this solution the 
points Xq and Xi are entirely disconnected: to a given configuration A°i{x) is associated a 
configuration of Euler angles at each point in space, {a{x), P{x),'j{x)), and to build Eq. (p^D, 
one simply reads off the values of these angles at Xq and Xi. 

But now it becomes clear that we have eliminated the previous problem at the cost 
of introducing another type of extraneous divergence, and this solution will also have to 
be discarded. The problem lies in the fact that the Euler angle variables appearing in the 
angular momentum operators, Eq. (|3^), originate in an integer-spin representation of SU{2), 
and thus take values in the SO {3) range < a < 27r, 0</?<7r, 0<7<27r (after all, they 
represent rotations of a tetrahedron in M^). It is in this 50(3) range that Eqs. ( p9| , |30D 
must be solved. However, by the time we solve them for j = 1/2 sources, Eqs. (|3^j35|) , we 
have had to resort to the fundamental representation of SU{2), which is a double- valued 
representation of S0{3), with range < a < Air, 0</?<7r, 0<7< 2tt. This means that when 
a (at either xq or xi) is rotated by 27r, which is no change at all in the SO (3) degrees of 
freedom, this solution changes sign. From the 5*0(3) viewpoint, the X>^/^ matrices undergo 
a discontinuous jump between a = 27r and a = 0. Of course, this has happened because local 
solutions such as Eq. ([3^ ) cannot be expected to automatically discern between SU{2) and 
50(3) degrees of freedom, a question which can only be decided at a global level. We are 
thus led to discard Eq. (|3lD as an appropriate solution to the Gauss law constraint with 
fundamental sources due to the double- valued and non-analytic nature of this solution. We 
also see, on the other hand, that the situation is entirely different with adjoint sources, where 
a solution would be Eq. ( ^41) with the substitution p(V2) _^ p(i) qj-, equivalently, 

^^^ = ut{x,)u%xoW'[g] . (36) 

It is trivial to see that this satisfies Eqs. (^,0) with A'^ in the adjoint and, in this case, 
\& would in fact be a single- valued function of the Euler angles in the 50(3) domain. We 
would furthermore be able to conclude that in this case no line appears connecting the points 
Xq and Xi, much like in the abelian case, and this seems to indicate there should not be a 
confining potential between these sources. At the same time, there are clear indications from 
lattice calculations that for a certain range of distances between the confinement and color 
screening scales, adjoint and higher representation Wilson loops do experience an area law 
behavior with a string tension proportional to the quadratic Casimir of the representation. 
While we are not questioning these results, what our results point out, however, is that one 
should exercise much caution in directly translating these Wilson loop scalings into the true 
potential experienced by static quarks. 



In attempting to find an appropriate solution to Eqs. (|2^ - |31|) with sources in the funda- 
mental, we are guided by the following prerequisites: i) we should not introduce all of the 
vector potential in a Wilson line prefactor, but rather only the "gauge part" of it, and ii) we 
should seek a functional which is an analytic function of the gauge variables in the 50(3) 
domain, so that the wavefunctional belongs to the correct Hilbert space. In separating the 
gauge part of the vector potential, the m" variables are again very useful. For, if we imagine 
that to each possible tetrahedron in space we associate a unique triplet of vectors f"(x), 
then the m" variables will in general be 50(3) rotations of these, 

<(x) = U^\x)v',{x) , (37) 



under gauge transformations. One must of course find a systematic way to assign a unique 
v^ to each tetrahedron, and once that is done, the purely gauge degrees of freedom become 
neatly separated in Eq. (0). Although it is not of central importance for our purposes here, 
we briefly describe how tetrahedra are associated to v^: given a tetrahedron, we align its 
longest edge with the z-axis with one vertex at the origin; we then align the longest edge 
connected to this first edge to lie on the left xz-plane (x < 0), also with one vertex at the 
origin. Then the last vertex will lie somewhere in IB?. The v? are then given by the three 
vectors describing the edges emanating from the origin, and some thought shows that the 
integration over gauge-inequivalent configurations corresponds to 

dri I dr2 / dr^i / d9i2 / dOi-^ / dcj) , (38) 

Jo Jo Jo Jo Jo 

where r^ = y/vfvf = ^/gu, 9ij = v'^Vj/riVj = gij/y/gligj] (no sum in i or j in these expres- 
sions), and is the azimuthal angle of fg (the vector not lying in the xz-plane). 

This separation of gauge- and gauge-invariant variations is sufficient to build the wave- 
functional we seek: it will be given by associating the [/"'', corresponding to a particular u^, 
to a pure gauge vector potential 

f/"^ -^ (A^Uix) = -^e'^'^U'"' {x)d,U'^'' (x) (r72)„, . (39) 

This vector potential is then used in constructing the Wilson line prefactor, leading to the 
candidate wavefunctional 

\ /a/3 

It is worthwhile to work out the gauge transformation properties of this wavefunctional in 
order to clarify the issues of analyticity and center symmetry. This is done by looking at 
the regularized form of its infinite product representation: 

N 

F-g= Ylil + tdx^Mx^,)) , (41) 

71=1 

where we understand the curve connecting xq and xi to be divided into a large number A^ 
of segments at the points {x„}, and Ai = A'}t°'/2. The first term in this product is 

(1 + tdx\Mxi)U = Sap - '-e'^'^U''' {xi)dU'''' {xi) f ^) , (42) 

with dU'^^{xi) = dx\diU"'^{xi). The standard commutation and anticommutation relations 
of Pauli matrices then lead to 

(1 + idx\A{x{)U = (yy) i^S"' + f/'^''(xi)f/''''(x2)) (43) 

to linear order, where X2 = Xi — dxi. One then uses the standard 2: 1 map between SU{2) 
and ^0(3), 
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''-'^T).,"-i-"'-[j)^,' '^^) 



- i5ai3'Si3a' " -^SafsSa'/B' j , (45) 



and the "spin-permutation" identity 

Z /a/3 V Z /a'/3' 

to finally write 

(1 + idx\A,{xr))^p = (f/(xi)f/t(x2))^^ [Ui}{U{xi)U\x2))] , (46) 

with tif being the trace in the fundamental. All A^ terms in the product can now be put 
together to obtain 

Kf = {U{x^)U\xo)U[ltTj{U{x^)UHx2))] . . . [^W(f/(x^)t/t(xo))] . (47) 

In the continuum limit one might be tempted to simply drop the infinite trace term since 
each trace involves two matrices which are very close together, and to linear order differ by 
an element of the algebra, which is traceless. This is almost correct, except for the fact that 
it endows the path ordered exponential with the 50(3) analyticity we have to require of 
state wavefunctionals. That this is so can be seen as follows: we would like our candidate 
wavefunctional, Eq. (|40|) , to be invariant under the S'0(3)-periodic change a ^ a + 27r, 
independently at Xq and Xi. This is equivalent to making a singular gauge transformation 
equal to the identity everywhere except for one of the points Xq or xi, where it is equal to 
a rotation of 27r around the z-axis. Then, either U{xi) or U{xo) will change sign, but since 
j?reg jg quadratic in the fundamental matrices at every point, it will remain invariant. By 
simply dropping the trace term, this periodicity would be broken. It is this analyticity in 
the 5*0(3) domain that we must preserve in taking the continuum limit of Eq. (^7]). Due to 
continuity, furthermore, the trace factor allows for a definition of the angles at Xi starting 
from a given convention for angles at Xq. Thus, if the angles at Xq are taken to be in the 
50(3) domain, then, in the continuum limit, the angles at Xi will be given by 

a'^(a;i) = q;(xo) + / dx'dialx), (48) 

and likewise for /3*^(a;i) and 'j'^{xi), where we write the superscript C over the angles simply 
to differentiate them from the angles obtained by simply reading them off at the point Xi, 
as in the previous solution (thus, angles at Xi have unlimited range). With this definition 
of angles at Xi, we are now ready to drop the trace factors and take the continuum limit for 
F: 

F^, = iU^ {x,)U\xo)U = 2)iV'H(«, A7)f)<^'^*((«,/?,7)o) . (49) 



This represents Eq. ( ^OD in a more explicit form. We note that it satisfies Gauss' law 
while at the same time being periodic in the 50(3) domain. This is achieved through a 
residual dependence on the line connecting the two charges. Yet, it is very similar to the 
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local solution Eq. (0), as it must be since both satisfy the same local gauge invariance 
constraint. Furthermore, it emulates the abelian case in the closest possible form, insofar as 
only gauge degrees of freedom are introduced to "communicate" between the two charges. 

It is worthwhile noting that if one were to form a fundamental Wilson loop in the fashion 
spelled out above, then it would transform with a Z2 factor, ±1, under gauge tranformations, 
with the specific sign being determined by the number of vortex-like singularities in the gauge 
transformation that pierced a surface bounded by the loop. This is most clearly illustrated by 
a gauge transformation with a single infinite vortex along the z-axis, which can be described 
as follows: 

o.{x, y, z) = tan~^ — , (50) 

X 

and (3{x,y,z) = 'j{x,y,z) = 0. Then, the Wilson loop will acquire a minus sign under this 
gauge transformation if and only if it encircles the z-axis an odd number of times. This 
shows how our analyticity requirement leads to a natural realization of 't Hooft's operator 
algebra |0 

A{C)B{C') = e2™/^5(C")A(C) , (51) 

where (in the notation of [@]) A{C) is a Wilson loop around C and B{C') is the operator 
which creates a vortex at C. For SU{2)/Z2 gauge group, the phase on the r.h.s. above is 
equal to the (±) Z2 factor described previously. If the gauge group were SU{2) rather than 
50(3) this Z2 factor would naturally be absent, and Wilson loops would be gauge invariant: 
in that case each point in space is uniquely associated to an SU{2) matrix and thus the 
matrices which transform the endpoints of a Wilson line become one and the same when the 
endpoints meet, and the transformation vanishes in the trace. 

It should also be noted that in order to achieve an expression for the wavefunctional 
of two sources which does not single out a preferred curve, one must naturally consider 
wavefunctionals which are linear combinations of the ones just described, going over different 
curves connecting the two charges, with some weighting factor w{C): 

^a, = T.^{C)F^,m- (52) 

c 

With such wavefunctionals at hand, one may of course then attempt to argue that confine- 
ment comes about due to random vortex fluctuations in the vacuum, as in the standard Z^r 
vortex condensation picture 0. At this point, however, we would rather refrain from such 
speculations. 

We finally observe that it is possible to exhibit the different wavefunctionals considered 
here by means of a description which is independent of the use of gauge covariant variables. 
If we imagine one were to select a family of gauge potentials. A", that crossed once and 
only once each gauge orbit, then any gauge potential A'} would be related to one of the A" 
through a gauge transformation (of course, trying to achieve such a family of A"; through 
a gauge condition inevitably leads to Gribov ambiguities, but as long as we do not choose 
gauge conditions we are free to assume the existence of the A" unambiguously). Then, the 
generic solution to the Gauss law constraint should, strictly speaking, be 

^a/3[A] = f/„^(xi)f/|„(xo)<l''""[i] , (53) 
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where f/^/j is the gauge transformation relating A'} and A?, and $™"[A] is to be determined 
by a dynamical calculation in the vacuum sector. The wavefunctional sometimes presented 
in the literature, gotten by acting with a fundamental Wilson line on the vacuum state is 

^«^[^]=P(e^/^-/2)^^ $.,,[!] (54) 

= U^m{x,)U;^{xo) [P (e^/^-/2)'"" ^^^[A]] . (55) 

It clearly suffers from the same problem as the incorrect wavefunctional in the abelian case, 
namely, the introduction of too many gauge quanta propagating between the two sources. 
Our solution, Eq. (|40|), would on the other hand, correspond to 

■^ap[A] = Uarnix,)U;^ixo) [6^^^' ^ ,^[A]] (56) 

= (f/(xi)f/(a;o)^)^^$vac[i] . (57) 

Of course, this description cannot be carried on much further, as we have no handle on how 
to describe the unique gauge slice A'}. 

IV. CONCLUSIONS 

We have investigated here the form state functionals should take in the presence of two 
static charges, both in abelian and in SU{2) gauge theory. Part of the motivation for such 
an investigation was the fact that the functionals which are suggested by the form Wilson 
loops take, are incorrect, and generate a confining static potential even where there is no 
confinement. In the abelian case, we have shown how the static Coulomb potential arises in a 
considerably more direct fashion through Schrodinger-picture methods than through Wilson 
loop expectation values. Another part of our motivation was to extend this treatment to the 
nonabelian theory, by exploring the gauge covariant description of the Hilbert space that we 
have proposed in previous work. In the nonabelian case, the use of gauge covariant variables 
has shown to be well-suited to the problem, as they allow for a direct solution of the Gauss 
law constraint before the dynamical problem is tackled. Their use, furthermore, together 
with an analyticity requirement, also led to a simple and explicit realization of the center 
symmetry of the theory. 

The extension to gauge group SU(3) should be essentially straightforward, whereas the 
actual dynamical problem of determining the static potential does present considerable tech- 
nical difficulties. We have begun such an investigation, and hope to report on it in the near 
future. 
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